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Abstract. We study modular determinantal differential equations of orders 2 
and 3. We show that the expansion of the analytic solution of a nondegenerate 
modular equation of type D3 over the rational numbers with respect to the 
natural parameter coincides, under certain assumptions, with the q-expansion 
of the newform of its spectral elliptic curve and therefore possesses a mul- 
tiplicativity property. We compute the complete list of D3 equations with 
this multiplicativity property and relate it to Zagier's list of nondegenerate 
modular D2 equations. 



1. Introduction 

Motivated by Apery's proof of irrationality of C(3), Don Zagier studies in [ZagierOT] 
the question of finding those triples of rational numbers {A, B, A) for which the se- 
quence obtained by the recursive formula 



(n + - (Air^ + An + A)u„ + Bn'^Un-i = 

starting with uq = I has all integer terms, i.e. u„ € Z. The generating function 
(l>o{t) = 1 + uit + U2t^ + . . . is the normalized analytic at t = solution of the 
differential equation £0o = with 

(1.1) C{t) = D'^ - t [AD^- + AD + \) + Bf{D + if , 

where we use the notation D = t-^ throughout the paper. We will refer later on 
to (jl.ip as the Beukers-Zagier differential operator since it appeared in the works of 
these two authors. A table of respective triples {A, B, A) is obtained in [ZagierOT] 
by searching in a large range of values. It appears that all degenerate cases in the 
table, i.e. those with either A^ = AB or B = come as members of infinite families 
of triples {A,B,X) with (l>o{t) E Z[[t]]. On the contrary, imposing the assumption 

(1.2) A"^ ^AB, S 7^ 

one arrives at 14 "sporadic" cases with no obvious pattern. Remarkably, in all 
those sporadic cases the corresponding differential equation can be parametrized 
by modular forms. Namely, one can find a modular function t(r) that vanishes at 
oo and a modular form /(r) of weight 1 such that (j)o{t{T)) = /(t) for all r in the 
upper half-plane with large enough imaginary part. These cases are listed in the 
table below. 



VASILY GOLYSHEV AND MASHA VLASENKO 



A 


B 


A 


t{r) 


f{r) 


"0 


Ul 


"2 


"3 


U4 


M5 


7 


-8 


2 


1^6" 
233s 


2^3" 
1263 


1 


2 


10 


56 


346 


2252 


9 


27 


3 


1^4^ 18^ 


2*^3"'^ 12"^ 


1 


3 


9 


21 


9 


-297 


2993363 


134303 


10 


9 


3 


283* 


2^3'^ 

1362 


1 


3 


15 


93 


639 


4653 


11 


-1 


3 


CO 

n=l 


/ 5\l/2 


1 


3 


19 


147 


1251 


11253 


12 


32 


4 


1*4=8* 
216 


210 
1*44 


1 


4 


20 


112 


676 


4304 


17 


72 


6 


1 = 3^4'*6=12^ 


2^^3=12= 


1 


6 


42 


312 


2394 


18756 


2l4 


164805 





-16 





2*8« 
41^ 


4IO 

2*84 


1 





4 





36 
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For each of the first 



The products in the forth and fifth columns stand for eta-products, e.g. 
ri{2T)^T](3T)B ■ notation throughout the paper. 

6 rows one should also consider {A, B, A) = {—A, B, —A) with the corresponding 
Un = (— l)"'Un, ^(t) = —t{T), f (t) = f (t) . For the last row there is also a triple 



(0, 16, 0) leading to u„ = (— 1) 



In total, the table gives us 14 triples {A, B, A) 



satisfying (|1.2p with (/)o(t) G '^[[t]]- Zagier conjectures that there are no more such 
cases, or if there are, they all will have modular parametrization. 

We observe in this paper that the differential operator p.ip satisfying the as- 
sumption (jl.2p is a specific form of the so called determinantal differential operator 
of order 2. The necessary definitions and properties will be recalled in Section [2j 
We deal with determinantal differential operators of order 2 in Sections |3] and |3] and 
recover Zagier's list above in a new context. Then we proceed to obtain an analog 
of this list for determinantal differential operators of order 3, namely, the complete 
list of D3's that satisfy a multiplicativity property that we discuss later. 

2. Determinantal differential equations 

Determinantal differential equations of order N were defined in |GS07] . A DN 
equation is obtained from an (A^ + 1) x {N + 1) matrix A — {aij)fj=Q that satisfies 

aij = , i — j > 1 
(2.1) a^i = 1, i-j = l 

aij — aN-j,N-i , * ~ J < 1 
The respective differential operator is then defined as 

where Sij is the Kronecker symbol and detright refers to the way of expanding the 
determinant of a matrix with non-commuting entries with respect to the rightmost 
column. 

The matrix A can be reconstructed from the coeflicients of the differential op- 
erator ( |GS07) . Corollary 3.3). Assume in addition that all eigenvalues of A are 
distinct. Then obviously A is diagonalizable. In fact, for a matrix satisfying (12. ip 
the two conditions are equivalent: A is diagonalizable if and only if all eigenvalues 
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of A are distinct. It follows immediately if one observes that A cannot have an 
eigenvector whose last component iz zero. According to Corollary 6.4 in |GS07j the 
singularities of the differential operator Ca,oo{z) are regular singular points located 
at oo and the eigenvalues Aq, . . . , Xn of A. Moreover, the differential equation has 
maximal unipotent monodromy at 2; = oo and the valuation of its analytic solution 
at z = cx) is equal to 1. This motivates the following notation. 

Definition 2.1. A differential operator of order N is of type DNoo.i if it equals 
i^A,oo{z) for some matrix A satisfying (|2.1[) . 

We denote the characteristic polynomial of A by F{z) — det(z — A) throughout 
the paper. It will be convenient to also use the variable t — j. Namely, consider 
the operator 

£aoW = {-If CA,oo(^)t. 

The respective differential equation has maximal unipotent monodromy aX t — Q 
and the valuation of its analytic solution at this point equals 0. 

Definition 2.2. A differential operator of order N is of type DNq^ if it equals 
^A,a{t) for some matrix A satisfying (|2.ip . 

By DN we mean either DNq^ or DN^o,! the case being clear from the context. 
Observe that the following operations with the defining matrices 

A^ A' = A + e 
A^A" = (A^-^+ifly) 
lead to the substitutions in the differential equations 

(2.2) Ca'Mz) = CA^oo{z-e), CA'fl{t) = £a,o (j^) (1 - 
and 

(2.3) La",M = A/:a,oo(^), LA"Ai) = ^Afi{\t) 
respectively. 



3. The Beukers-Zagier equation as a D2 equation 
Let us consider D2 equations in detail. According to our definitions one has 




where 

F{z) = Aet[z - A) = + a2Z^ + aiz + ao 
a2 = -an - 2aoo 



ai = 2aooaii + apo - 2aoi 
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A Z)2oo,i differential operator is then any operator of the form 

with a cubic monic polynomial with distinct roots F{z) — + a2z'^ + aiz + ao. 
One can recover the matrix A from ai and /3 via 

aoo = /? 

Oil — —2/3 — a2 

3 2 1 

ooi = "2/^ - /3q!2 - -ai 

0-02 = - ^^^2 - Pai - ao 
The generic equation of type -D2o,o is then 



tJKdtJ \ \tJ \tJJdt \t 

d\2 



with 

= t + a2t^ + ait^ + aaf^ 

H{t) ^ - 13 + ait + 2aQt^ 

With the notation D — t-^ we can further rewrite it as 

(1 + a2t + aif + a^t^){D'^ - D) + {I + 2a2t + Sait^ + 4aoi^)£' 

(3.1) - /3i + ait^ + 2ao<^ 

= + t{a2D^ + aal? - /?) + ait2(i:) + l)^ + aot^{D + 1){D + 2) 

Notice that putting ao = we obtain precisely operator (jl.ip with A = —a2, 
B = ai and \ — (3. 

4. Modular equations D2 

Recall that a D2 differential equation depends on 4 parameters (a2j cki, aoj /?)• 
It determines a local system of rank 2 over the base 

P^(C) \ {00, the roots of z^ + a2Z^ + aiz + ao} ■ 

Consider also the basis in the space of solutions of D2 near t = which consists on 
normalized analytic and logarithmic solutions: 

</.o(<) = l + /3<+(-ia2/3 + i/32-iai)i2 + ... 

</.i(t) = \ogt^a(t) + (-a2-2/3)i+(ia2 + ia2/3-^/32- iai)i2 + ... 

Definition 4.1. M^e say that an equation D2 with parameters (a2, ai, ao, /3) G 
is modular if the analytic continuation of 

1 ^i{t) 

T — 

27ri (/'o(i) 

gives uniformization of the base by the upper halfplane with the group of deck trans- 
formations being a congruence subgroup o/SL(2,Z) and the function t i— >■ (/>o(t(r)) 
is a modular form of weight 1. 
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In this case ^(t) is a modular function whose q-expansion can be written exphc- 
itly. Indeed, inverting the series 

q = exp(^) = t + {-a^ - 2P)f + [al + ^a^/? + f Z?^ - \a,)t^ + ■ • ■ 
one gets 

9 19 1 

t = q + {a2 + 2p)q^ + [al + -a2p + — + -ai)g^ + . . • 

Further, substituting this expansion into 4'o{t) one obtains 

/ = Mt{r)) = 1 + /3<Z + (i«2/3 + J/32 - + . . . 

This must be a modular form of weight 1. 
Put Q = q^ and consider the series 

00 

(4.1) ^3 ^1 + + + aot^ (j)o{tf = ^ c„Q" 

n=l 

whose coefficients c„ = Cn{d, (3) can be determined explicitly as follows. One writes 

^ - MlW)^ = tHl + (-^"^-/3)^+(^2 + 5«./3 + ^/32-l«,)t^ + ...) 
and inverts this series in order to get 

= g (l + (ia2+/3)Q2 + (^a2 + 1^2/3 + y/3' + iai)Q4 + ...) 

which can be substituted into the left-hand side of (|4.1|) . We have 

ci = 1 
C2 = 
C3 = ^2 + 3^ 

(4.2) C4 = 

2 25 ^ 75 . 1 



It is not hard to see that all even coefficients in fact vanish. 

Theorem 4.2. Assume one is given a nondegenerate modular D2 with parame- 
ters a2,OLi,oio, P G Q • Consider the weight 2 modular form X^J^Li Cn?" with c„ 
determined from the expansion (|4.1I) . // in addition it is a newform then 

00 

L{s) = y^ 

n=l 

is the L-function of the elliptic curve 

(4.3) y'^ — + 02^2 + aiz + ao ■ 

By being a newform we mean that the modular form belongs to the subspace on 
newforms of certain level. We do not require it to be a Hecke eigenform a priori; 
rather, the Hecke-eigen property is a consequence of the theorem. In particular, 
we have the following 

Corollary 4.3. // X^riLi '^"9" ^■^ newform, then its coefficients (|4.2p are multi- 
plicative, i.e. 

(4.4) Cmnia,l3) = Cm(d, /3) • c„(a,/3) 
as soon as m and n are coprime. 
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We will solve equations (|4.4p with respect to the parameters a2,ai,aoj/3 later 
in this section. It appears that modulo a certain transformation which preserves 
both the sequence {c„;n > 1} and the L-function of (j4.3p there are finitely many 
cases. 

The proof of Theorem 14.21 will rely on the following result. 

Theorem (Atkin & Swinnerton-Dyer congruences, Theorem 4 in |ASD71) ) Let 

p ^ 2, 3, and let = + Bz + C he an elliptic curve over Zp with good reduction. 
Choose a local parameter at so that z = + Xl^-i '^n?" o.'^d y — ^^"^ + . . . 
are the respective expansions , and write 



Idz 



n=l ^ 

If B,C,dmCn are p-adic integers, then 

(4.5) Cnp - OpCn + PC2^ = mod 

where 

^ /m^ + Bm + C\ 

«p = - Li )■ 

Notice that this theorem can be applied to an elliptic curve defined over Q with 
good reduction at p as soon as the coefficients B, C, dn and c„ do not contain p in 
their denominators. Moreover, Op = p + 1 — ^E(¥p) is then the p-th coefficient of 
the L-function of this elliptic curve. 



Proof of Theorem \4.2\ Let a„, n > 1 be the coefficients of the L-function L{s) 
J2n Tfi of the elliptic curve (|4.3p . One can check that 

^[it)Mt) - Mm'oit) = ,J ,2^ ,3^ = ^"'^(t 

t(l + a2t + ait'' + aot-^) \t 

hence 

= i^F(i)*(*;(<)*o(i)-*i(0*iW)<if = \t-'F(^y*dt = -i^ 

where we substitute z — 1/t, y'^ — F{z). This is a holomorphic differential on the 
curve (|4.3p . and since Q ^ t^ for small t we conclude that Q is a local parameter 
on the curve near the origin. Moreover, z ~ and y ~ and therefore 

the theorem of Atkin and Swinnerton-Dyer would be applicable for every prime p 
not dividing the conductor of the curve as soon as all c„ and d„ defined from the 
expansion z — j — + dnQ" do not contain p in their denominators. 

First we show that this is indeed the case for all but finitely many primes p using 
the assumption of modularity. 

We have Q{t{T)) = . Looking at (|4.ip we see that X^^i "^n?" ^he q- 
expansion of the modular form t^\/l + a2t + ait'^ + aot^f^ (2t) of weight 2. It 

follows that possibly after multiplication by an integer all c„ become integers si- 
multaneously. The same holds for dn since z{t) is a modular function. Therefore 
for all but finitely many prime numbers p we have congruences (|4.5p . Another 
consequence of modularity of Cn?" is that 

(4.6) c„ — o{n) , n — > oo . 
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Our next step is to show that (|4.5p together with (|4.6I) imply that c„ = a„ for all 
n not divisible by a finite set of primes. Since a„ = o(n2+^) for any e > and 
c„ = o(rt) there is a nmnber iV such that 

1 



<2 



1 

<2 



for all n > A^. Obviously we can assume that (|4.5p is true for all p > increasing 
N if necessary. From (j4.5p with ti = 1 we get 

Cp = Up mod p . 

Since for all p > TV also |cp|, |ap| < | we conclude that Cp = Op. Suppose p > N 
and we have proved that Cpm = Opm for all m < M. Since 



ttpA/ + l — CLpCLpM -\~pCLpM-l — 
2pCp 



CpM+1 — apCpM +pCpM-i = mod p^^^"'^ 



we conclude that 

CpM+i = CpM+i mod p'^'^'^^ . 

Therefore CpM+i = UpM+i because p^'^+^ > N again. It follows now by induction 
that Cpm = Cpm for all m > and p > N. Our next step is to show that Cn = a„ for 
all n not divisible by finite number of primes p < N. Let n be such a number and 
suppose that for every proper divisor n'\n we have already proved that a„' = c„'. 
By (j4.5p with a prime divisor p\n and n' — ^ instead of n we have 

Cn — cipCn' + pc^ — c„ — a„ = mod p™ 

where m — ordp(n). Since this is true for every prime dividing n we conclude that 
Cn = a„ mod n, and therefore c„ = a„ again by our estimate. 

Consider both newforms J2n Cn?" and J2n '^nl^ on the intersection of the cor- 
responding congruence subgroups which is again a congruence subgroup. Since 
Cn = o-n for all n not divisible by primes form a certain finite set, it follows that 
both forms have the same eigenvalues for infinitely many Hecke operators. There- 
fore by multiplicity one theorem these forms are just equal, and our theorem is 
proved. □ 

Now we can substitute the polynomials (|4.2p into the multiplicativity relations (|4.4p 
and solve the resulting equations. We do not expect finitely many solutions because 
the shifts (j2.2p preserve modularity. Under the shift A' = A+e the parameters 
become 

(a2, a'l, ttp' /'O — («2 — 3e, 3e^ — 2a2e -I- ai, — -I- a2e^ — cue + Q;o, /3 + e) • 

Zagier's choice ag = is not natural from this point of view because one can make 
— Q only if F{z) has a rational root. It is more natural to choose the equation 

with 13 — Q as & unique representative of the orbit of the shifts. 

Solving the first few relations with Grobner bases ( we used computer algebra 

system [MAGMA, ) 

C6 = C2 • C3 , Cio = C2 • C5 , C12 = C4 • C3 , C14 = C2 • C7 , 
Cl5 = C3 • C5 , C18 = C2 ■ Cg , C21 = C3 • C7 , C22 = C2 ■ Cn 

we obtain 8 points (a2,ai,ao) plus two one-parametric families (0,0, ag) and 
(0, ai,0). In order to show that there are actually finitely many cases in these 
families we used more relations by considering about 200 further coefficients. The 
results are given in the table below. 
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a2 


ai 


ao 


I 


n 

yj 


n 

W 


-1 








-4 


-80 


-192 


4 


-80 


192 


-2 


-40 


-75 


2 


-40 


75 


-1 


-24 


-36 


1 


-24 


36 





16C4 











27C6 



Fjz) 

z^{z + l) 

z2(z- 1) 
z- 12)(z + 4)2 
z + 12)(z-4)2 



z + 3)(z^-5z-25) 
z-3)(22 + 5z-25) 
z-6)(z + 2)(z + 3) 
z + 6){z - 2){z - 3) 
z(z2 - I6C4) 
- 27C6 



-1,0,0),(2,1,1) 
1,0,0), (-2,1,-1) 
-32, 256, -12), (16, 0,4) 
32,256, 12), (-16, 0,-4) 



11,-1,3) 
-11,-1,-3) 

-17, 72, -6), (7, -8, 2), (10, 9, 3) 
17,72,6), (-7, -8, -2), (-10,9, -3) 
0,±16,0), (12,32,4), (-12,32,-4) 
9, 27, 3), (-9, 27, -3) 



The first three columns contain all solutions of a few first multiplicativity equa- 
tions. The fourth column shows the roots of the respective polynomial F{z) = 
z^ + OL^z^ -f (X\z -\- QfQ. Polynomials in the first 4 rows appear to have multiple roots 
meaning that the respective differential operator is degenerate. In the last column 
we shift the differential operator by various roots of F(z) in order to obtain opera- 
tors with ao = 0, the respective values of the parameters (A, i?. A) — (— a2,cn,/3) 
being listed. The last 6 rows give us D2 equations, and shifting by various ratio- 
nal roots we obtain precisely Zagier's table. One can easily check at this point 
that in each case the statement of Theorem 14 . 2 1 holds . and therefore the respective 
(a2,Q!i,ao) indeed solve all multiplicativity equations. The triples correspond- 
ing to the degenerate differential equations from the first four rows can be found 
in |Zagier07| as ^^1, #3, #19 and #11. On the other hand, the degenerate triples 
#14, #20 and #25 are also modular but do not appear on our list. 



5. Differential equations of type D3 

Our goal in this section is to write the generic form of a D3 equation by making 
exactly the same steps as in Section [3] but now with = 3. We get 



F 



det ri 



ght 



l{z-a,,)± -aoi(f)' -ao2{±f 

-1 {z-a,,)£ -aoi(^) 

\ -1 {z-aoo)£j 



with 



F{z) = det(z-^) ^ z^ + 



asz 



a2Z 



aiz 



as 
a2 



-2aii — 2aoo 



4aooaii -|- ang -|- a,. - 2aoi - ai2 



^dz> 



ai = -2ao2 - 2aiiaQQ + 2aoo(aoi - Oil + 012) -f 2aiiaoi 
"0 = -ao3 + 2aooao2 + {aj^ - 012)090 - 2aiiaoiaoo + 
G{z) = z^ + piz + pQ 
I3i -2aoo 

/^o = 2aooaii - aj^ - 2aoi + ai2 

Recall that a 15300,1 differential equation is called non-degenerate whenever the 
roots of F{z) are distinct. Notice that our order 3 differential operator is the 
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symmetric square of the order 2 operator 



We also compute the generic -DSo.o- We have 



= t 
with 



Hit){ 



dt 
di' 



-F' 



d\2 



dt) 
1 



t - 



F" 



-r 



dt) 



a2t'^ + aii^ + a^f 



Hit) . 

U{t) = i^G(i) + 3t"*F(i) ~ ^^^'(t) = Saot^ + ^ait^ + {a2 + + I3i 



Finally, this differential operator can be written as 



(5.1) 



+ t{D + ^){a3{D^ + D) + /3i) 
+ t\D + l){a2iD + lf + /3o) 



a 



,tHD + 2){D + -){D + l) 



+ aat\D + 3){D + 2)(D + l) 

6. NONDEGENERATE MODULAR EQUATIONS D3 



In this section we will prove the analog of Theorem l4.2l for D3 equations. In order 
to state it we first associate to such an equation an appropriate elliptic curve. Recall 
that F{z) = det(z — A) has distinct roots, so the discriminant of F is nonzero. For 
a D3 equation we have F{z) = + asz^ + a2Z^ + aiz + ao- Consider the 



z** + asz^ + a2z'^ + a\z 



curve 
(6.1) 

Put it into the Weierstrass form. 

Definition 6.1. The spectral elliptic curve of a D3 equation is 



1 2 

y = x" + (aia3--a2 



4ao I X 



/ 2 1 2,8 
+ [oioa^ - -aia2a3 + — a2 - -aoa2 + c 

It is indeed an elliptic curve because the discriminant of the cubic polynomial in 
the right-hand side is equal to the discriminant of F{z) (as a function of «»), and 
therefore the right-hand side has 3 distinct roots. 

Lemma 6.2. Curves (|6.ip and (|6.2p are birational over the splitting field of the 
polynomial F[z). Moreover, the holomorphic differential ^ on the spectral elliptic 
curve transforms into on (|6.ip under this birational equivalence. 

Proof. Let F(zo) = 0. Then 

— (z — zo)* 4- az{z — zo)^ -I- a2(z — zq)^ + ai(z — zq) 



with Ui = ^F('')(zo), and 
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Hence the variables 

6.3 X = h — , y 

satisfy = x'^ + Bx + C with 

Q!2 1 2 

S = aiQ!3 — ^ aiQ!3 - -q;2 - 4ao , 

C = <3? + ~ ^<5ia2a3 = "oaa - ^aia2a3 + ^"2 ~ ^aoa2 + Oi\ . 

The equahty of the differentials follows immediately. □ 

Consider the normalized analytic and logarithmic solutions of D3 near t = 0: 
Mt) - 1 - \ht + (^"3/31 + ^/S? - - i/?o)i' + . . . 
^^■^^ Mt) = \ogtMt) + i'(t) - logt0o(<) + (-ia3 + ^/3i)t + ..• 

Again, the 6 parameters naturally split into two groups. Parameters as, a2, ai , ao 
determine the base 

P"'"(C) \ {00, the roots of + asz'^ + a2Z^ + a\z + ao} 

and also the spectral elliptic curve (j6.2p . Of course the solutions (|6.4p and the 
respective local system of rank 3 depend also on the remaining parameters ji\ and 
/3o. 

Definition 6.3. We say that an equation D3 with 03, a2, ai, aoj j^ii /3o G Q *s 
modular if the analytic continuation of 

gives uniformization of the base by the upper halfplane with the group of deck trans- 
formations being a congruence subgroup o/SL(2,Z) and the function t i— >■ (pQitir)) 
is a modular form of weight 2. 

Consider the power series 

We can invert this expansion in order to write t as a power series in 5, and for 
accessory values of parameters this must be then the g-expansion of a modular 
function. Analogously, 4>o{i) written as a power series in q will be the q-expansion 
of a modular form of weight 2 since the local system is of rank 3. We will be 
specifically interested in the coefficients of the modular form 

oc 

(6.6) t^t) = ^c„g" 

n=l 

which can be computed explicitly as polynomials in the initial parameters 
ci = 1 

C2 = ^^3 - Pi 

^^■''^ 3 2 27 „ 57^2 1 3 ^ 
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Theorem 6.4. Assume we are given a nondegenerate modular D3. If the modular 
form (j6.6p 0/ weight 2 is a newform then 



00 




n=l 



is the L-function of the spectral elliptic curve (|6.2p . 



Proof. We will use the same method as for D2. Let a„,n > 1 be the coefficients 
of the L-function L(s) = ^" " ^^'^ spectral elliptic curve (|6.2p . One can 
check that 

tcj)o[t) — = t(po{t)d\ogq = t — dt 



(6.8) 



Fi'iyK-'dt = 

t J w 



where we substitute z — 1/t, — F{z). Now using the birational transforma- 
tion (j6.3p we see that t (hence also q) is a local parameter on the spectral elliptic 
curve near the point P = ^^Q!2, &ij where ~ jjF^^^zq) and zq is a chosen root of 
F{z) = as in the proof of Lemma 16.21 Since our differential equation is modular, 
composition of modular uniformization with the birational transform gives a map 
from a modular curve to the spectral elliptic curve. The preimages of both points 

±P = (^^a2,iaij are cusps because t is zero there and this is a cuspidal value as 

we know. Therefore by the Manin-Drinfeld theorem their difference P— {—P) = IP 
and hence also P is a point of finite order. Let us write the spectral curve (|6.2p 
as y'^ = + Bx + C. The differential — ^ in (|6.8p transforms to ^ according to 
Lemma [121 One can find an isogenous curve = i'^ + Bx + C where P is mapped 
to the origin and the differential is mapped to —j^- Then q is a local parame- 
ter on the latter curve near the origin, and the L-function of this curve is again 
L(s) — On n^^ since isogenous curves have equal L-functions. One can check 

that since the expansion — {^^=1 Cnq"'^ ^ starts with ci = 1 then x ^ q~^. 

All the coefficients c„ in (j6.6p and dn in the expansion x = q^^ + do 
not contain p in denominators for all but finitely many primes p because these are 
g-expansions of a modular form of weight 2 and modular function respectively. The 
rest of the proof goes exactly like in Theorem 14.21 Namely, one has Atkin and 
Swinnerton-Dyer congruences for all but finitely many primes and together with 
c„ — 0(71), which is another consequence of modularity, this implies that a„ — Cn 
for all n non divisible by a finite set of primes. The normalized newforms a,nq^ 
and "Ylin Cn are then in one Hecke-eigenspace and therefore are equal by multiplicity 
one theorem for the space of newforms. □ 

By this theorem, if X^i^i Cn?" is a newform then the coefficients c„ are the 
coefficients of the L-function of an elliptic curve, and we have the following conse- 
quence. 

Corollary 6.5. //X^i^i Cn?" is a newform, its coefficients (|6.7p are multiplicative, 
i.e. 

(6.9) Cmn{a.,fi) = Cm(d, /3) • c„(a,;3) 

as soon as m and n are coprime. 

Let us solve equations (|6.9p for the parameters {a, (3). 
Lemma 6.6. A D3 equation has the following properties with respect to the .shift 
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(i) the parameters change according to the rule 



(a3'"2>"'l>"o) 

— (as — 4e, a2 — Seas + 6e^, ai — 2ea2 + Se^as — 4e^, ao — eai + e^a2 — e^^aa + e^) 
(/3;,/3^) = (/3i-2e,/3o-/3ie + e2) 



(ii) the coefficients of the spectral elliptic curve (j6.2p do not change 

(iii) all On in (|6.7p do nof change 

(iv) modular differential equations transform into modular ones (assuming e € 

Q) 



Proof. Indeed, since for A' = A + e one has Ca',oo{z) — Ca,oo{z — e) the formulas 
in (i) follow from 



z"* + a^jZ"^ + a2Z^ + a'lZ + a'^ 

= (z - e)* + a^{z - e)^ + a2(-z - ef' + ai(z - e) + ao, 
z^ + l3[z + P'o = (z-e)2+/?i(z-e) + /?o. 



Part (ii) follows by a tedious computation. Next, from CA'.o[t) ~ 'Ca.o j-rjj j (1~ 

et) we conclude that Frobenius bases must simply transform as (t)'i{t) — (1 — 

£t)4>i{^YZ^^ for i = 0, 1,2. Therefore the uniformization maps (|6.5I) differ by the 

transformation i j-rjj of the base, from where (iii) and (iv) follow immediately. 

□ 



It follows from Lemma 16.61 that it suffices to solve equations (j6.9p for a single 
representative of every orbit of shifts. There are two natural choices, as = and 
/?! = 0. We will use the latter one. It appears that the system of equations (|6.9p 
has finitely many solutions with /3i = 0. We will give the complete list later, but 
first we list the rational solutions that correspond to nondegenerate DS, i.e. such 
that the roots of F{z) are all distinct. There are exactly 18 of them. We list the 
ai 's In the table below. These determine the spectral elliptic curve (j6.2p which we 
denote by £. Then we give its j-invariant, its level N and the newform g[T) of level 
N whose Mellin transform is the L- function of S. We give the value of /3o in the 
last column. 
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as 


a2 


ai 


ao 




N{£) 


9£(T) 


Po 





-44 





-16 


20720464 
15625 


20 


2^ • 10^ 


-4 





44 





-16 


80 


4® • 20^/2^ • 8^ • 10^ • 40^ 


4 






-28 
28 






-128 
-128 


207646 
6561 


24 

48 


2 • 4 • 6 • 12 
4'* • 12'*/2 • 6 • 8 • 24 


-4 
4 





-40 





144 


35152 
9 


24 


2 • 4 • 6 • 12 


-8 





40 





144 




48 


4* • 12^/2 • 6 • 8 • 24 
4* • 12^/2 • 6 • 8 • 24 


8 





-68 





1152 


'^nfi'^fii 71 

OUDODl 1 LO^ 

9 


48 


-28 





68 





1152 


24 


2 • 4 6 12 


28 





-48 





512 


287496 


32 


42 . 82 


-16 


n 


4:0 


n 

u 








42 . c2 


1 fi 











-256 


1728 


32 


42 .§2 














256 




64 


8^/42 . 








-36 





432 


54000 


144 


I212/6* • 24* 


-12 





36 





432 




36 


6* 


12 


-4 



-88 



-300 
-108 


-304 



122023936 
161051 




11 
27 


l2 . Il2 

3^ 9^ 


-8 




-2 


-43 


-156 


-216 


4733169839 
3515625 


15 


1 • 3 • 5 15 


-5 


-2 


-59 


-136 


-80 


4956477625 
941192 


14 


1 • 2 • 7 • 14 


-5 



In the proof of the Theorem 16.41 we have constructed the point P = (^■^0:2, (Xi^ 
of finite order on the spectral elliptic curve. The order of P is 4 in the first 14 cases 
and for the last 4 rows it is 5,3,8 and 3 respectively. 

7. All solutions of the multiplicativity equations for D3 

The goal of this section is to list all the solutions to the multiplicativity equa- 
tions. We have obtained them via Grobner bases with the aid of the computer 
algebra system [MAGMA] , doing computations over several finite fields and lifting 
solutions afterwards. Apart from the non-degenerate cases which we listed in the 
previous section, there are solutions defined over number fields and also solutions 
with degenerate F{z). In order to list them in an efficient way we consider the 

twists A I— >■ A" — ^A^^^+^flyJ which lead to the simple variable change in the 

differential equation CA".o{t) = CA.oi^t). 

Lemma 7.1. A D3 equation has the following properties with respect to the twist 
CA"fl{t) = CAfl{\t): 

(i) the parameters transform according to the rule 

a" = (Aa3,A2a2,AV,A'^ao), /3" = [XPuX^M 

(ii) the spectral curve transforms via 

^ + X^Bx + X^C 

(iii) the coefficients Cn{d,P) transform via c'( = X^~^Cn 



(iv) for the function T{t) 



1 0i(*) 

2TTi 0o(t)' 



ne has 
log A 



27ri 



The proof is straightforward. The way the solutions transform under the twists 
described in (iii) and (iv) of Lemma 17.11 shows that only finitely many twists are 
possible for a given D3 that preserve the multiplicativity property, and the only 
twists possible are those by roots of unity. We give the complete list of solutions 
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to (|6.9p below. We list only one representative in every family of twists and give 
all possible twists in the right-most column. We start with nondegenerate cases. 
These have already been given in the previous section up to the twists by roots of 
unity. 



as 


a2 


ai 


ao 


m 


N{£) 




/3o 


twists 





-44 





-16 


20720464 

2» 
6561 


20 




2^ . 10^ 


-4 


A2 = ±1 





-28 





-128 


24 


2 


4 • 6 • 12 


-4 


A2 ^±1 





-40 





144 


35152 
30656^17154 


24 


2 


4 • 6 • 12 


-8 


A2 = ±1 





68 





1152 


24 


2 


4 • 6 • 12 


28 


A2 ^±1 





48 





512 


287496 


32 




42.32 


16 


Ai«^ = l 











-256 


1728 


32 




42 .§2 





Ai'^ = l 





36 





432 


54000 


36 




6^ 


12 


A3'^ = l 


-4 


-88 


-300 


-304 


122023936 
161051 


11 




l2 . Il2 


-8 


A = ±1 








-108 








27 









Ai8 = l 


-2 
-2 


-43 
-59 


-156 
-136 


-216 
-80 


4733169839 
3515625 

4956477625 
941192 


15 
14 


1 
1 


3 5 15 
• 2 • 7- 14 


-5 
-5 


A = ±1 
A = ±1 



In addition, there is one more nondegenerate solution over Q(\/5) 



as a2 ai ao 


m 


/3o 


twists 


22-30V5 1000-440\/5 




18- 10\/5 


A2 = ±1 



which appears to give the same modular form i(r)(/)o(T) = 2^ ■ 10^ as in the first 
row of the above table. 

The solutions of (j6.9l) with degenerate polynomial F{z) = 2^ + a^z^ + 022;^ + 
aiz + ao are given in the table below. Remarkably, some of them are still modular. 



as 


a2 


ai 


ao 


/3o 


t{T)Mr) 


twists 

















q 







4 








-4 


q/{l-q^) 


A2 = ±1 





-8 





16 


-8 


g/(l-<z2) 


A2 = ±1 





128 





4096 


64 


2* • 8'*/44 
48/2^ 


A2 = ±1 





64 











A2 = ±1 


4 














qlil^qf 


A = ±1 


2 


1 








-1 


- q) 

ql{l + q + q^) 


A = ±1 


-2 


-3 








3 


A = ±1 


-6 


-135 


-540 


-648 


-9 


1^ 93/32 
24.64/12.32 

33 . 6^/1 • 2 

5Vl 
2^0 -8^/1^ -46 


A = ±l 


4 


-80 


192 





-16 


A = ±l 


2 


9 


-216 


432 


-9 


A = ±1 


2 


-55 


-100 


1000 


-25 


A = ±l 


8 


-176 


768 


-1024 


-16 


A = ±l 



8. From D2's to D3's 

Let 4'o{t) = 1-1- u\t + U2t'^ + ■ ■ ■ be the solution of the differential equation p.ip 
with ao = 0. As we already mentioned, this is exactly the class of D2's considered 
by Zagier in |Zagier07] where his parameters A, B, A are our — a2, ai,/3 correspond- 
ingly. One then has 

{n + l)^u„+i + (a2rt^ + a2n — (i)un + ai7i^u„_i — , 
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and we observe that the sequence it^ — (^,")wri satisfies 

(n + l)^w„+i + 2(2n + l)(a2n^ + - /3)m„ + 4ai(2n + l){2n - l)nun-i = . 
In other words, (/)o(i) = J2^=o (^")""^" ^ solution of 

+ t{D + i) {ia2{D^ + D) ~ 4/3) 
+ i^{D + l)(l6ai(i:» + if - 4ai) 
which turns out to be a D3 equation with parameters 

(a^,4,a;,a;,) = (4a2, 16ai, 0, 0) , (/?;, /3^,) = (-4/3, -4ai) . 
This equation is degenerate as its symbol has double roots. Passing to (t)'^{t) = 
(^r)'""^^"' ■^hich is a solution of 

+t^{D + l){Aa2{D + 1)2 - 16/3 - Aa^) + 16ait''(i:' + i){D + 2)(L> + 1) , 

we thus come to a D3 equation whose parameters are 

(4',a^',a'/,0 = (0,4a2,0,16ai), (/3;', /3^') = (0, -16/3 - 4^2) • 

This D3 is nondegenerate if the initial D2 was nondegenerate. Indeed, the symbol 
F"{z) = + Aa2z'^ + 16ai has 4 distinct roots if and only if ai 7^ and 7^ 4ai, 
which are exactly the conditions for the roots of the symbol F{z) = + a2z'^ + aiz 
of the initial D2 to be all distinct. Therefore we have a map from nondegenerate 
D2's with ao = to nondegenerate D3's with 013 = ai = /3i = given by 

d={a2,ai,0) Q? = (0, 4q!2, 0, 16ai) 

(8.1) . ^ ^ 

/3=(/3) /3 = (0,-16/3 -4a2) 

This map is obviously bijective. The analytic solution at f = transforms as 



n=0 



E 

n=Q 



Unt 



We find that under the map (|8.ip Zagier's triples all go into D3's from our list 
of accessory equations and, moreover, they exhaust all D3's in our list with = 
ai =0. There 14 cases on both lists, below we show half of them as we already 
did in Introduction. Namely, for every D2 listed in the table with the parameters 
(q;2, Q;i, ao) there is also (— 02, ai, — ao). We denote the spectral elliptic curves of 
the D2 and D3 differential equations by £2 and £3 respectively. 







D2 










£13 






Q!2 


ai 


P 




N{£2) 


Q!2 


ao 


/3o 


jX^a) 




-7 


-8 


2 


1556068 
81 


24 


-28 


-128 


-4 


207646 
6561 


24 


-9 


27 


3 





144 


-36 


432 


-12 


54000 


144 


-10 


9 


3 


1556068 
81 


24 


-40 


144 


-8 


35152 
9 


24 


-11 


-1 


3 


488095744 
125 


20 


-44 


-16 


-4 


20720464 
15625 


20 


-12 


32 


4 


1728 


32 


-48 


512 


-16 


287496 


32 


-17 


72 


6 


1556068 
81 


48 


-68 


1152 


-28 


3065617154 
9 


48 





-16 





1728 


32 





-256 





1728 


32 
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The spectral elliptic curves appear to be non-isomorphic in general but their levels 
coincide. In fact, the respective spectral curves are isogenous, and therefore their 
L-functions are equal. Indeed, the spectral curve of D2 

f 2 : = + a2z'^ + aiz 

is isogenous to the elliptic curve 

f 2 : 2/^ = - 2q;2Z^ + (q!2 - 4q;i)z , 

the isogeny of degree 2 given by 

'y2 y[ai-z^)^ 



{z,y) ^ 



z^ ' 



(see Example 4.5 in |Silvernian09] ). This latter curve is in turn isomorphic to the 
spectral elliptic curve of the respective D3 



2 3 / 16 2 X .128 3 512 

£3 : r = z-" + (-ya2 - 64ai)z + (^"2 - —otxa^) 



with the map 



References 



[Zagier07] Zagier, Don, Integral solutions of Apery-like recurrence equations, Harnad, John (ed.) 
et al., Groups and symmetries. From Neolithic Scots to John McKay. Selected papers of the 
conference, Montreal, Canada, April 27-29, 2007. Providence, RI: American Mathematical 
Society (AMS) (ISBN 978-0-8218-4481-6/pbk). CRM Proceedings and Lecture Notes 47, 349- 
366 (2009). 

[GS07] V. Golyshev, J. Stienstra, Fuchsian equations of type DN, Commun. Number Theory 

Phys. 1 (2007), no. 2, 323-346 
[ASD71] Atkin, A.O.L.; Swinnerton-Dyer, H.P.F., Modular forms on noncongruence subgroups, 

Combinatorics, Proc. Sympos. Pure Math. 19, 1-25 (1971). 
[Silverman09] Silverman, Joseph H., The arithmetic of elliptic curves., 2nd ed. Graduate Texts 

in Mathematics 106. New York, NY; Springer, 513 p. (2009) 
[MAGMA] W. Bosma, J. Cannon, and C. Playoust, The Magma algebra system. I. The user 

language, J. Symbolic Comput., 24 (1997), 235-265 

Algebra and Number Theory Sector, Institute for Information Transmission Prob- 
lems, B. Karetny 19, Moscow 127994, Russia 
E-mail address: golyshevamccme.ru 



School of Mathematics, Trinity College, Dublin 2, Ireland 
E-mail address: masha.vlasenkoagmail.com 



